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PREFACE 


This  report  was  prepared  at  the  Georgia  Tech  Engineering  Experiment 
Station  under  Contract  No.  DAAB07-74-C-0272.  The  work  covered  by  this 
report  was  performed  in  the  Applied  Engineering  Laboratory  under  the  super 
vision  of  Dr.  H.  A.  Ecker  and  Mr.  J.  L.  Eaves,  Director  of  the  Applied 
Engineering  Laboratory  and  Chief  of  the  Radar  Technology  Division,  respec- 
tively. The  progress  reported  herein  was  performed  during  the  eighth 
quarter  of  a program  to  develop  an  environment  and  radar  operation  simula- 
tor (EROS)  to  be  used  in  testing  radar  receivers  and  components. 


This  project  is  being  monitored  by  Mr.  Reinhard  G.  Olesch  and  Mr 
Otto  E.  Rittenbach  of  the  U.  S.  Army  Electronics  Command,  and  their 
helpful  suggestions  are  acknowledged. 
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1.  INTRODUCTION 


This  report  covers  work  performed  during  the  period  1 April  1976 
through  30  June  1976,  the  eighth  quarter  of  a 30-month  program  to  design 
and  build  an  Environment  and  Radar  Operation  Simulation  (EROS).  The  effort 
has  been  devoted  primarily  to  continuing  the  implementation  of  EROS 
hardware  and  software. 

A large  part  of  the  EROS  hardware  is  completed  through  board  design, 
wiring,  and  component  assembly.  The  following  major  units  have,  in  addition, 
been  completed  through  unit  test  and  debug: 


Computer  Interface 

This  unit  moves  data  from  the  PDP-11  unibus  to  various  locations 
in  the  EROS  digital  hardware.  It  provides  protocol  for  conversing  with 
PDP-11  I/O  instructions  and  contains  a small  amount  of  buffer  storage. 


IK  by  16-bit  Memory  Module 

The  IK  by  16-bit  memory  module  is  used  in  the  clutter  filter,  the 
target  processor,  and  in  the  video  signal  generator.  The  memory  module 
includes  drivers  for  addressing  and  enable  lines. 
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Clutter  Timing  and  Control 

This  unit  generates  the  control  sequence  which  defines  the  opera- 
tions performed  by  the  clutter  hardware.  It  synchronizes  the  clutter  and 
target  processors,  and  it  controls  the  filters,  azimuth  weighting,  and 
azimuth  Integration. 


Clock  Generator 

The  clutter  generation  hardware  receives  4 KHz  and  10  MHz  clock 
signals  from  the  analog  hardware.  The  4 KHz  clock  is  derived  from  the 
radar's  P.R.F.  The  10  MHz  clock  is  derived  from  a crystal  oscillator. 
The  clock  generator  unit  provides  drivers  of  these  signals,  and  it 
divides  the  10  MHz  clock  by  four  yielding  four  phases. 
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Hardware  Clutter  Test 


This  unit  provides  facilities  for  troubleshooting  and  spot  checking 
the  operation  of  the  clutter  generation  hardware.  It  is  capable  of 
forcing  the  clutter  hardware  to  stop  when  any  specified  memory  location 
is  addressed,  and  it  permits  single-step  operation  from  the  stopped  state. 
It  also  produces  test  data  for  exercising  the  clutter  processor  independent 
of  the  PDP-11  and  of  other  EROS  hardware. 


A significant  portion  of  the  simulation  preparation  and  real-time 
software  has  been  completed.  The  following  routines  have  been  tested 
and  debugged: 


Spectral  Clutter  Library  Maintenance 


This  collection  of  programs  maintains  a library  of  clutter  spectral 
information.  Each  spectral  record  in  the  library  is  identified  for  reference 
purposes  by  an  integer  called  a "spectral  number".  These  programs  diagnose 
spectral  library  updates  to  detect  user-generated  errors,  derive  filter 
statistics,  and  maintain  a reference  file  in  spectral-number  order.  The 
derivation  of  filter  statistics  is  described  in  Section  2 of  this  report. 


Complete  Clutter  Libracv  Maintenance 


This  collection  of  programs  maintains  a library  of  complete  clutter 
!i:ormation.  Each  clutter  record  in  the  library  contains  spectral  and  radar 
cross-section  information  and  is  identified  by  a two-character  name.  These 
programs  diagnose  clutter  updates  to  detect  user-generated  errors,  derive 
fill  or  parameters,  and  maintain  a reference  file  of  complete  clutter 
.at.)  motion  in  identifier  order. 


Range- Dependent  Reference  Data 


Each  simulation  run  uses  a table  of  64  weights  representing  the 
range-dependent  radar  signal  attenuations  (e.g.  R”^  power  loss).  These 
64  weights  are  selected  from  a table  of  197  entries  corresponding  to 
60  m.,  75. m . . . 3000  m.  The  particular  entries  selected  depend  upon 
the  band  of  ranges  being  simulated.  This  program  diagnoses  the  input 
data  to  detect  user-generated  errors  and  translates  the  data  into  12-bit 
numbers  for  real-time  simulations. 


Ta rget  Recording 


These  programs  collect  (complex)  samples  of  recorded  target  data 
from  a dual-channel  A/D  converter  and  store  the  digitized  data  on  PDP-11 
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disk  for  future  reference.  After  the  recording  is  finished  the  signal  is 
reread  from  the  disk  to  calculate  its  mean  squared  amplitude  for  the 
purpose  of  radar-cross-section  scaling.  Since  the  samples  are  read  at 
a 4 KHz  rate,  one  of  the  programs  operates  with  real-time  constraints. 


Real-Time  Executive 

This  program  keeps  track  of  time  and  antenna  position.  It  invokes 
target  processing  of  the  illuminated  target  (if  any),  starts  clutter 
input/output  when  appropriate,  and  permits  the  execution  of  any  non-time- 
critical  functions  (e.g.  CRT  display) 


Real-Time  Target  Processor 

This  program  maintains  target  azimuth,  obtains  antenna  weights 
from  a core  table,  and  transmits  this  data  with  target  samples  and  range 
information  to  the  EROS  digital  hardware.  The  real-time  requirements  for 
this  program  are  critical,  and  the  code  is  optimized  for  speed.  The 
disk  reading  portion  of  this  program  (to  fetch  target  data)  has  not  yet 
been  written,  but  the  program  has  been  tested  and  debugged  with  data 
initially  generated  in  core. 
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2.  CLUTTER  FILTER  TRUNCATION  STATISTICS 


; 
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The  physical  realization  of  the  clutter-synthesis  filter  introduces 
errors  due  to  the  finite  length  of  digital  number  representation.  These 
errors  are  called  truncation  errors,  and  they  occur  specifically  when 
two  16-bit  numbers  are  multiplied  to  yield  a 16-bit  product.  This  section 
contains  a brief  review  of  the  source  of  truncation  errors  followed  by 
a description  of  the  algorithms  that  compute  their  effect  on  the  statistics 
of  the  filter  outputs. 


2.1  Source  of  Truncation  Errors 

Recall‘d  that  the  complex  digital  filter  which  synthesizes  clutter 
has  two  delay  elements,  whose  outputs  at  time  nT  (after  the  nth  iteration) 
are  denoted  by  the  complex  quantities  w^(nT)  and  W2(nT0.  The  calculations 
performed  by  the  filter  are  described  by  the  equations 


wx  ^ (n  + 1)T 
w2  ( (n  + 1)T 


) 

) 


w (nT)  + ot2  w2(nT)  + x(nT) 
w1  (nT) 


v(nT) 


6 w (nT)  + K 


(1) 


where 


x(nT)  is  the  sequence  of  complex  random  Inputs; 

and  a 2 are  (real)  parameters  that  control  the  clutter  spectral 
statistics; 

S is  the  (real)  parameter  that  controls  the  radar  cross  section 
(RCS)  of  the  moving  portion  of  the  clutter  in  the 
simulation  cell; 

K is  the  (complex)  parameter  that  controls  the  RCS  of  tl§.  immobile 
portion  of  the  clutter  in  the  simulation  cell; 

v(nT)  is  the  sequence  representing  the  sampled  complex  envelope 
of  clutter  return. 


This  formulation  is  extracted  from  Section  2 of  the  Sixth  EROS 
Quarterly  Report. 


I 


i 


The  filter  is  being  implemented  in  the  form  of  two  identical  modules; 
one  generates  real  components  of  v(nT)  for  the  544  cells  illuminated  within 
an  antenna  beam,  and  the  other  generates  imaginary  components.  For  reasons 
discussed  in  Section  2.5  of  the  5th  EROS  Quarterly  Report,  the  product  of 
d]W^  is  actually  formed  as  the  sum  (a  - l)w^  + wj..  The  numbers  in  the 
digital  filter  modules  are  represented  as  16-bit  2's  complement  binary 
fractions  with  the  assumed  binary  point  between  the  sign  bit  and  15 
fraction  bits.  Thus  during  every  computation  cycle  the  real  and  imaginary 
components  of  the  two  products  (a^  - l)w^  and  a2  w 2 are  truncated  from  30 
fraction  bits  to  15  fraction  bits. 

The  errors  introduced  by  the  truncation  are  expressed  as  follows. 

For  any  real  number  a,  let  [a]  denote  the  largest  integer  that  does  not 
exceed  a.  For  any  complex  number  a + jb,  let  [a  + ib]  denote  [a]  + j [b] . 
Then  the  truncated  products  are  [2^  (a^  - l)w^]2_1'  and  [2  5 a 2 w2]2“  5, 
and  the  truncation  errors  are 


el  = t2l5(ai  " 1)wl]  2"15  ‘ (al  ~ 1)wl 

and  (2) 

e2  « [215  «2  w2]2-15  - a2  w2 

The  effects  of  truncation  errors  can  be  analyzed  by  treating  and  &2 
as  additional  sources  of  random  input,  because  Equation  (1)  is  still 
valid  if  x(nT)  is  considered  to  be  the  sum 


x(nT)  * e1  (nT)  + e2  (nT)  + r(nT) 


(3) 


where  r(nT)  is  the  output  produced  by  the  pseudo-random  number  generator. 


2.2  Relevance  of  Truncation  Statistics 

2 

For  any  random  variable  u,  let  y{u}  and  s {u}  denote  its  mean  and 
variance  respectively.  If  u is  complex  with  real  part  uR  and  imaginary 
part  Uj,  then 


p{u>  ■ + Jvfuj) 

and 

s2{u}  = s2{ur}  + S2{Uj} 


(4) 


1See  Section  2.3  of  the  6th  EROS  Quarterly  Report. 
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Using  this  notation,  the  third  equation  in  (1)  can  be  restated  in  the  form 
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v(nT)  = Btw^nT)  - y {wj.}]  + By^}  + K (5) 

The  first  term  in  (5)  has  zero  mean  and  represents  the  radar  return  from 
the  moving  scatterers  in  the  clutter  cell.  The  second  term  and  third 
term  in  Equation  (5)  are  constant,  and  their  sum  represents  the  radar 
return  from  immobile  scatterers  in  the  cluttjr  cell.  The  RCS  of  the 
moving  scatterers  is,  therefore  proportional  to  the  mean  squared  ampli- 
tude of  B[w^  (nT)  - y{wj_}],  which  equals  B2s^{w^}.  Similarly  the  RCS  of 
the  immobile  scatterers  is  proportional  to  | By{w^}  + k|2.  Thus  in  order 
to  calculate  6 and  K from  the  moving  and  immobile  RCS  values,  it  is 
necessary  to  estimate  y{w^}  and  s^{w^}.  These  quantities,  in  turn,  can 
be  computed  from  the  following  formulas. 


F = 


G = 


1 ~ al  " “2 


1 - a- 


(1  + a2) [ (1  - a2)2  - Oj2] 


(6) 


F[(l  + j)  y{p}  + y {e1>  + y (e2>] 


s2{w1> 


F2 (~2yy  {p}  - 2[y{p}]2)  + G2 (2y2  + s2{e1>  + s2{e2>  + 2yy{p}) 


Mt.  ! 

I j 

J 


The  parameter  y is  the  magnitude  of  the  pseudo  random  input  and  is  selected 
as  the  largest  number  from  the  set  {2  , 2-2,  . . . , 2”  5}  for  which  the 
filter  is  safe  from  overflow  (i.e.  for  which  |Re{w.  (nT) } | < 1 or 
|lm  (w^(nT)}|  < 1).  The  quantity  y{p}  equals  y(231  - l)-*.  From 

these  relationships  it  follows  that  the  calculation  of  3 and  K from  the  two 
clutter  RCS  values  requires  that  y,  y{e^}  , s^{e^},y{e2],  and  s ie2)  be 
evaluated. 


It  is  clear  from  symmetry  that  for  each  of  the  random  variables 
x(nT),  e-^  (nT),  (nT),  and  w^  (nT),  any  property  which  holds  for  the  real 


^Refer  to  Section  3 of  the  7th  EROS  Quarterly  Report. 
2Refer  to  Section  2.3  of  the  6th  EROS  Quarterly  Report. 
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part  of  the  random  variable  also  holds  for  the  imaginary  part.  Let  the 
subscript  R denote  "the  real  part  of"  and  the  subscript  I denote  "the 
imaginary  part  of".  Then  for  any  of  the  above  random  variables  u 


\ 

I 


p{u}  = (1  + j)  * (1  + j)  ytuj.} 

S2{u>  = 2s2  {Uj^}  = 2s2{ui>  . 


(7) 


2.3  Evaluation  of  y 

To  calculate  y the  following  inequality  is  employed: * 

00 

| w1R(nT)  | < | xR  | I h(nT)  I ’ (8) 

n=l 

where  | xR  | is  the  maximum  of  | x^  (nT)  | and  h(nT)  is  the  filter's 
impulse  response.  Let  e^R  and  e2R  denote  maximum  values  of  | e^R  (nT)  | 
and  | e2R  (nT)  | respectively.  Then  clearly  if 

A A °° 

(Y  + 81R  + e2R>  1 h<nT)  1 < 1 (9) 

n=l 

the  filter  will  not  overflow  with  pseudo-random  input  of  magnitude  y. 

The  values  of  e^R  and  e2R  are  easily  calculated  from 

e1R  = 2‘15  (1  - 2_mi) 

(10) 

^2R  = 2“15  (1  - 2_m2) 

where  m^  and  m2  are  the  smallest  respective  integers  such  that  2™1 (a^  - 1) 


1c.  f.  Oppenheim,  A.  V.,  and  Schafer,  R.  W.,  Digital  Signal  Processinj 
Prentice- Hall,  Englewood  Cliffs,  N.  J.,  1975,  p.  427. 


Is  somewhat 


and  2m^x2  are  integers.  The  evaluation  of  ^ ' | h(nT)  | 

more  complicated.  The  denominator  of  the  z-Sransform  of  h(nT)  is 
2 

z - a^z  - and  if  the  polynomial  either  has  two  identical  roots  or 

00 

has  0 for  one  of  its  roots,  then  a closed  form  expression  for  ^ ^ | h(nT) 
can  be  derived.  In  the  two  identical  roots  case  n"^ 


£ 

n=l 


h(nT) 


1 - 


when  + 4a2  - 0 (11) 


In  the  0-root  case 


00 

Z|  h (nT)  | = when  a„  = 0 

1 - I a,  I 2 


(12) 


Although  no  closed-form  expression  has  been  found  for  the  general 
case,  a numerical  evaluation  of  the  sum  is  straight  forward.  To  approximate 
the  sum  the  following  expression  for  a bound  on  the  remainder  has  been 
found  useful 


£ 


h(nT)  | 


m 


£ 


n=l 


2 a 


m + 1 


h(nT) 


max 


(1  - a ) 
max 


(c^  + 4a2  T*  0) 


+ 4a„ 


(13) 

where  a is  the  magnitude  of  the  larger  of  the  two  roots: 
max 

amax  " 2 max  { | °1  ± V “l2  + 4a2  | } (14) 

A computer  program  to  calculate  y using  these  formulas  has  been  written  and 
debugged . 
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2.4  Evaluation  of  Truncation~Error  Means  and  Variances 


An  iterative  approximation  method  is  being  employed  to  estimate  the 
means  and  variances  of  e^  and  • Observe  that  Equation  (2)  expresses 
the  random  variables  e^  and  as  functions  of  w^  and  W2*  where  - 1 and 
Oj  are  constants.  We  define 


eQ  (w)  = [ 215  aw  J 2~15  - 


so  that  e.  * e 1 (w, ) and  e0  - e (w_).  Since  the  only  difference 

1 - .1  1 t Oj  L 

between  w.  and  w2  is  a delay  of  one  step,  the  probability  distributions 
of  w1R  and  w2R  are  identical.  Therefore,  if  their  density  function  p(w1R) 

were  known,  then  the  means  and  variances  of  the  error  variables  could  be 
computed  from 


{ ea  ("lR>  } ' y:  % <W1R>  P <"lR>  U‘> 

W1R 

{ea  *W1R*  } ■ 'y  v ea  (wlR5  P (wlRJ  " [y{*a  (wlR^  (17) 


1 


It  is  necessary  to  treat  P(w1R)  as  a discrete  probability  distribution 
function  (and  to  use  sums  instead  of  integrals),  because  for  most  of  the 
values  of  a,  ea  (w.R)  is  not  even  approximately  a continuous  function  of 
w^R.  The  discrete  set  of  values  assumed  by  w^R  is  {+i2“15:  i * 0,  . . . 

215  - 1). 

Unfortunately  P(w1R)  is  not  known.  However,  computer  simulations 
have  empirically  shown  that  the  distribution  of  w.R  is  approximately 
Gaussian^-.  Therefore,  we  can  assume  the  approximation 


P(i2"15) 


2tt  s{w1R} 


(i2~ 15  - U(w._})2 


2s^wiR} 


^efer  to  Section  2.4  of  the  Fifth  EROS  Quarterly  Report. 
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and  iteratively  estimate  p{w  } and  s (w  }.  The  iterative  algorithm 
is  more  specifically  the  following: 

1.  Initialize  y{e1R}  and  y{e2R}  to  ^ eiR  and  2 e2R  resPectlvely* 
(Refer  to  Equation  (10).  The  set  of  possible  values  assumed  by  e^  is 

{ 0,  1 • 2 ^ mi,  . . . , (2mi  - l)2  ^ mil.  The  largest  value  in 

A 

this  set  is  e._,  and  the  average  value — assuming  equal  likelihood  for 

A*'  ^ A A 

each  possibility  — is  -r  e1D.  The  same  discussion  applies  to  e_.) 

2 ^2  ^ 

Initialize  s {e^R}  and  s {e2R)  to  zero* 

2 

2.  Calculate  bfw^R)  and  s {w^R}  using  (6). 

3.  Recalculate  p{e1R>  and  s2{e1R}  using  (15),  (16),  (17)  and  (18). 

4.  A heuristic  criterion  for  convergence  is  applied:  another 

iteration  starting  at  Step  2 is  performed  if  the  newly  computed  values 

of  p{e1D}  and  s^{e.  } are  not  sufficiently  close  to  the  previous  estimates. 
1R  IK 

The  computation  of  eQ(w^R)  is  somewhat  easier  than  that  suggested 
by  Equation  (15).  The  following  algorithm  is  employed  in  the  program 
which  calculates  means  and  variances. 


1. 

Compute  r^. 

residue  (modulo  2^)  ^ of  2^a. 

2. 

Compute  r2. 

the  residue  (modulo  2^)  of  w^R. 

3. 

Compute  r^. 

the  residue  (modulo  2^)  of  r^r2. 

4*  ea(wlRJ  2 r3 


^For  any  two  Integers  m and  x,  the  residue  (modulo  m)  of  x is  the 
Integer  r such  that  x ■ qm  + r,  where  q is  an  Integer  and  0 £ r < m. 


3.  MAXIMUM  RADAR  CROSS  SECTION  FOR  CLUTTER 


3.1  Introduction 


The  radar  cross  section  corresponding  to  an  EROS  signal  is  pro- 
portional to  its  average  squared  amplitude.  More  specifically 


where 


a is  the  RCS  of  the  simulated  scatterers, 

o 

K is  the  proportionality  constant  (20,000  m for  AN/PPS-15 
simulations),  and 

~ 2 

A is  the  average  squared  amplitude. 


The  components  of  RCS  associated  with  a given  clutter  cell  are 
denoted  by  oR  and  os  for  the  moving  and  immobile  scatterers  respectively. 
Since  their  respective  mean  squared  amplitudes  are  62s2{w^}  and 

| K + 6y{w^}  | (refer  to  the  previous  section  of  this  report), 

- Kc  b2“2<»1)  - 2Ko  b2°2<”ir> 


os  - Ko  |K  + Mw^l2  - K0  |k  + Mw1r}(1  + j)  r 


3.2  Parameter  Constraints 

The  magnitude  of  the  parameter  8 has  an  upper  bound  v,  which  is 
determined  by  the  position  of  the  assumed  binary  point  in  the  digital 
hardware.  If  there  are  n bits  between  the  sign  and  the  assumed  binary 
point,  then  v - 2n.  In  the  current  implementation  n - 0;  therefore, 
v ■ 1.  Since  oR  is  proportional  to  82  as  expressed  by  (20), 


Htefer  to  the  Section  3 of  the  Seventh  EROS  Quarterly  Report. 
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°R  * 2*o  v2  s2  '»1R> 


(21) 


The  real  and  imaginary  components  of  K are  also  limited  by  the 
positions  of  the  assumed  binary  point.  This  binary  point  position  has 
been  selected  to  be  consistent  with  that  of  B.  Thus 


H 

i 

' 


< v 


(22) 


Kx  | < v 


These  inequalities  lead  to  a mutual  constraint  between  oR  and  Og. 

The  computer  program  CCDIAB  maintains  the  library  of  reference 
data  for  different  types  of  clutter.  The  user  supplies  spectral  and  RCS 
information,  and  CCDIAB  stores  this  data  in  a file  for  subsequent  retrieval. 
CCDIAB  calculates  B from  oR  using  the  first  formula  in  (20). 


2K„  {.lR> 


(23) 


However  CCDIAB  does  not  have  sufficient  information  to  calculate  K in 
similar  fashion  because  the  phase  of  K is  not  determined  until  scenario 
compilation  time.  The  clutter  scenario  compilation  program,  which  utilizes 
the  reference  file  produced  by  CCDIAB,  randomly  assigns  the  phase  of 
S ■ K + Bv(w^}  from  a uniform  distribution  between  0#  and  360°,  and 


assigns  /og/Ka 


to  the  magnitude  of  S in  accordance  with  (20). 
with  S defined,  K can  be  calculated 


Then 


K - S - By{Wl} 


(24) 


In  order  to  satisfy  (22) 


|SR  - Mw1R}|  < v 


and 


(25) 


Sj  - Mv1R}|  < V 
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Utilizing  the  fact  that  | y{w  } | < 1 and  that  g < v,  it  can  be 

shown  that  the  inequalities  (25)  are  true  for  all  complex  numbers  S of 

magnitude  J Og/Ko  if  and  only  if 

aS  < Ko  (v  " 6 MW1R>D2  (26) 


Combining  (23)  and  (26) 


(27) 


Inequality  (27)  expresses  a mutual  constraint  limiting  the  magnitude  of 
the  radar  cross  sections  aR  and  Og. 

The  following  is  a summary  of  the  computations  and  tests  performed 
by  CCDIAB,  which  ensure  that  the  pertinent  constraints  are  met. 

1.  The  radar  cross  section  a D of  the  moving  scatterers  is  compared 

2 2 K 

to  the  maximum  value  2K  v s {w  } implied  by  (21).  If  o_  is  too  large, 

O J.K  K 

it  is  reduced  to  the  maximum  value. 

2.  6 is  calculated  from  aR  using  (23). 

3.  The  radar  cross  section  Og  of  the  immobile  scatterers  is  compared 

to  the  maximum  value  K (v  - 6|y{w  }|)^  implied  by  (26).  If  o_  is  too 

O IK  J> 

large  it  is  reduced  to  the  maximum  value. 


3.3  Dependence  of  Maximum  RCS  on  Spectral  Parameters 

The  radar  cross  section  of  the  immobile  clutter  scatterers  is,  for 
all  practical  purposes,  non-constrained.  Values  of  |v(w^R}|  tend  to  be  very 
small  (much  less  than  .5).  .Therefore,  the  maximum  valueof  Og  that  satisfies 
(26)  is  greater  than  5000  tn  for  worst-case  g (6  ■ 1),  and  5000  m2  is  only 
about  6 dB  below  the  largest  EROS  RCS  (20,000  m2). 


The  radar  cross  section  of  the  moving  scatterers  is,  on  the 
other  hand,  sensitive  to  the  spectral  parameters.  Maximum  values  of  aR, 
for  typical  cutoff  frequencies  (f_)*  calculated  from  (21),  are  tabulated 
in  Table  I. 

•j  i 
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TABLE  I. 


Maximum  Radar  Cross  Section 


(•2) 


Cutoff  Frequency  f o max 

K 


20 

211 

15 

124 

10 

100 

5 

12 

2 

Values  of  s {w  } were  estimated  by  the  computer  program  which  uses  the 
algorithm  described  In  Section  2.  The  spectral  parameters  and  In 
these  calculations  were  derived  from  two-pole  analog-filters  whose  s-plane 
poles  subtend  + 150°  with  the  p sitive  real  axis  and  have  magnitudes 
appropriate  for  the  given  cutoff  frequency. 


2 

Table  I Illustrates  the  fact  that  s {w^}  decreases  quickly  with 
fc.  This  is  not  surprising  since  the  dominant  term  in  the  Formula  (6) 

for  s2{w^}  is  2G2  y2,  and 


g2  y2  » Zh2  (nT) 

(E  | h(nT)  |)2 


(28) 


o 

(the  numerator  equals  G by  definition).  As  fc  decreases,  the  number  of 
non-negllglble  items  in  h(nT)  increases,  and  the  ratio  in  (28)  is  therefore 
reduced . 
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4.  NEXT  QUARTER  PLANS 


The  results  presented  in  Section  2 and  3 represent  the  completion 
of  the  analysis  to  support  EROS  design.  During  the  next  quarter  the  effort 
will  be  devoted  primarily  to  implementation,  and  design  additions  or 
modifications  will  be  performed  only  if  omissions  or  errors  are  identified. 
It  is  anticipated  that  all  of  the  EROS  hardware  will  have  been  completed 
by  the  end  of  next  quarter  through  unit  assembly  and  test.  As  time 
permits  part  of  the  hardware  integration  may  be  performed.  The  plans 
for  software  development  next  quarter  include  clutter  scenario  definition 
via  CRT,  target  scenario  definition  via  CRT,  clutter  scenario  compilation, 
and  real-time  clutter  processing. 
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